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Entanglement is a striking feature of quantum mechanics and an essential ingredient in most 
applications in quantum information. Typically, coupling of a system to an environment inhibits 
entanglement, particularly in macroscopic systems. Here we report on an experiment, where dis- 
sipation continuously generates entanglement between two macroscopic objects. This is achieved 
by engineering the dissipation using laser- and magnetic fields, and leads to robust event-ready 
entanglement maintained for 0.04s at room temperature. Our system consists of two ensembles 
containing about 10^^ atoms and separated by 0.5m coupled to the environment composed of the 
vacuum modes of the electromagnetic field. By combining the dissipative mechanism with a contin- 
uous measurement, steady state entanglement is continuously generated and observed for up to an 
hour. 

PACS numbers: 



To date, experiments investigating quantum superpo- 
sitions and entanglement are hampered by decoherence. 
Its effects have been studied in several systems [T] . How- 
ever, it was recognized [1] that the engineered interaction 
with a reservoir can drive the system into a desired steady 
state. In particular, dissipation common for two systems 
can drive them into an entangled state [3]. The idea of 
using and engineering dissipation rather than relying on 
coherent evolutions only, represents a paradigm shift with 
potentially significant practical advantages. Contrary to 
other methods, entanglement generation by dissipation 
does not require the preparation of a system in a partic- 
ular input state and exists, in principle, for an arbitrary 
long time, which is expected to play an important role 
in quantum information protocols [IHZ]- These features 
make dissipative methods inherently stable against weak 
random perturbations, with the dissipative dynamics sta- 
bilizing the entanglement. 

We report on the first demonstration of purely dissi- 
pative entanglement generation [5]. In contrast to pre- 
vious approaches [SHU], entanglement is obtained with- 
out using measurements on the quantum state of the en- 
vironment (i.e. the light field). The dissipation-based 
method implemented here is deterministic and uncondi- 
tional and therefore fundamentally different from stan- 
dard approaches such as the QND-based method [S] or 
the DLCZ protocol [3], which yield a separable state if 
the emitted photons are not detected. Furthermore, we 
report the creation of a steady state atomic entanglement 
by combining the dissipative mechanism proposed in |12j 
with continuous measurements. The generated entan- 
glement is of the EPR type, which plays a central role 



in continuous variable quantum information processing 
[3 US] , quantum sensing [T3] and metrology [TTl [THl |TB] . 

Fig. la presents the principles of engineered dissipation 
in our system consisting of two -'^^^Cs ensembles, interact- 
ing with a y-polarized laser field at wl. A pair of two- 
level systems is encoded in the 6S1/2 ground state sub- 
levels It)/ = |4,4)j, |;), = |4,3)/ and |t)// = |4, -3)//, 
|i)„ = |4, -4)//. Operators J±jj with J- = ^^I, |t).(i| 
describe collective spin flips, where N is the number of 
atoms. The atoms are placed in a magnetic field in the x- 
direction and the collective operators Jy — \/2 ( J+ -I- J~) 
and = iy/2{J+ - J") are defined in the frame rotat- 
ing at the Larmor frequency fl. The two ensembles are 
initialized by optical pumping along the x-axis in the ex- 
treme states mp = A and mp = —4 respectively, corre- 
sponding to (Jx) = {Jx,i) = -{J'x.ii) ~ 4A^ (see Fig. 1). 
Within the Holstein-Primakoff approximation, we intro- 
duce the canonical variables Xj/u = Jyj/n / \/\{Jx)\ 
and Pi/ii = ±J^j/ulyJ\{Jx)\ EI- The EPR entangle- 
ment condition [51 [T7] for such ensembles is given by ^ = 
Sj/ (2|(J,)|) = var(X/ - Xii)/2 + var(P/ + P„)/2 < 1, 
where Sj = var(Jj,j — Jyji) + var(Jzj — Jz.ii)- The 
entangling mechanism is due to the coupling to the a:- 
polarized vacuum modes in the propagation direction z 
of the laser field (Fig. 1), which are shared by both ensem- 
bles and provide the desired common environment. Spin 
flip processes in the two samples accompanied by forward 
scattering result in indistinguishable photons leading to 
quantum interference and entanglement of the ensembles. 
These spin flips and the corresponding photon scattering 
(see level schemes in Fig. 1) are described by the inter- 



2 




(b) Experimental setup: 

^Sn-detectors 

□ □ »^ 

pump-beams ^2,cos 




. 852nm 
pumping ; 
scheme: , /;-t''893nm 



FIG. 1: a) Collective dissipation modes and atomic levels: 
two spatially separated atomic ensembles interact with the 
environment composed of the vacuum modes of the electro- 
magnetic field. The coupling is driven by the j/-polarized laser 
beam. The engineered collective dissipation is due to photons 
scattered in the forward ^-direction. Internal level scheme of 
the atoms: the effective two-level systems |t) and |4,) are two 
magnetic sublevels Zeeman-shifted by a magnetic field ap- 
plied in the a;-direction, which defines the quantization axis. 
Atoms in the two ensembles are initialized in opposite spin 
states. The laser beam off-resonantly couples these levels to 
the excited states and to the electromagnetic vacuum modes. 
Due to the Zeeman shift of the ground state levels, photons 
are emitted into the upper and lower sidebands (shown in 
blue and red color) leading to collective spin flips J±. b) Ge- 
ometry of the experiment. The S2 detector signal processed 
by the lock-in amplifier (LI) is used to determine the atomic 
quantum spin components Jy,^ as described in the text. The 
optical pumping scheme is also shown. 



action Hamiltonian of the type 

H « /a.,, (^«k + ^^flk) + dk i^Bai + i?tak) 

where the integrals cover narrow bandwidths centered 
around the lower and upper sideband at =F f2 re- 
spectively and with the non-local spin operators A — 
fiJf — yJ^, B = /iJj^ — yJ^ . The fact that the electro- 
magnetic modes form a continuum is crucial for the 
entanglement to be created without measurements |12j . 
As emission into the forward direction is collectively en- 
hanced for a large optical depth d the forward scat- 
tered modes can successfully compete with spontaneous 
emission modes in directions other than z which leads to 
decoherence of the atomic state. Note that the genera- 
tion of entanglement cannot be explained by the inter- 
action of photons emitted by the first ensemble with the 
second one, which is negligible in our parameter regime. 
The nonlocal dissipative atomic dynamics obtained after 
tracing over the photonic modes is governed by the mas- 
ter equation [Hj: 

j-p = d ^ {ApA"^ - A^Ap + SpSt - B^Bp + H.C) + 
'CnoisoP, where p is the atomic density operator, and F is 



the single atom radiative decay. The Lindblad terms in 
parentheses, which would usually describe regular spon- 
taneous emission, drive the system into an EPR state 
with ^ — {^i — < 1 [TJ|, due to the special nonlo- 
cal construction of A and B. £noiso describes undesired 
processes such as single atom spontaneous emission, col- 
lisions, etc. The experiments are performed using two 
dilute ^'^■^Cs gas samples in 2.2cm cubic cells separated 
by 0.5m described elsewhere [B]. A bias magnetic field of 
0.9G leads to a Zeeman splitting of = 27r • 322kHz 
(see Fig. 1). The anti-relaxation coating of the cell 
walls and careful magnetic shielding [Tl] provide the non- 
radiative decoherence time for populations and coher- 
ences of Ti w 130ms and T2 ~ 40ms. The two ensembles 
are initialized in the states |4, ±4) with orientation up 
to P = 0.998(3) by applying a pump laser polarizing the 
F = A manifold and a laser repumping atoms from F = S 
to F = 4 for 10 to 50ms (Fig. lb). The driving laser is 
blue detuned by 850MHz from the F = A ^ F = 5 tran- 
sition of the D2 line corresponding to (/i — y)'^ — 0.16. 
The laser power influences both the collective and the 
single atom dissipation processes and has been optimized 
within a range of 5 — 15mW. The nonlocal atomic state 
variance ^ = J^j / {2\{Jx)\) is inferred, and the entangle- 
ment condition ^ < 1 is verified by a local polarization 
measurement on the light transmitted through the two 
ensembles (Fig. lb). We use the same laser to create and 
to verify the entanglement which significantly simplifies 
the experiment. In the period i < T, up to a variable 
time T (see the pulse sequence in Fig. 2b) the laser serves 
only as the driving source for dissipation. The results of 
the measurements on the transmitted light are not used, 
which is equivalent to tracing out the light field. Begin- 
ning at t = T, the temporal mode of the transmitted 
light is used for the determination of the atomic state at 
time T using the established method [SI [TTl [TSl [H] of 
linear mapping of the atomic state onto light (atomic to- 
mography via quantum polarization spectroscopy). The 
particular linear mapping used here has been utilized in 
several other contexts [TU [191 EQ] and is described by 
the input-output relations for atomic and light operators 
before and after the interaction: 
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and similarly for Xi — Xji Pi + Pu and yc ~> Us- Here 
— [\~ e^'^'^''^)/{fj, — v)"^ and 7^ oc (^ — i/)^ Ja;$, where 
$ is the flux of photons in the drive field and T is the 
interaction time. The light operators are given by the 
cos(r2t) and sin(Ot) components of the Stokes operator 
^2 weighted with an exponentially falling (rising) mode 
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out(in) 



rl QOut(il 
Jo "^2 



{t) cos(m)e=F'r=*dt 



-/+ Af±yS^JO '-^2 

(analogously for sine modes). Just as the master equa- 
tion does, the input-output relations predict an entan- 
gled atomic state with variance S,/(t) — )• {p — y)'^, for 
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FIG. 2: Entanglement generated by dissipation (a-c) and steady state entanglement (d). a) Time evolution of Ej(f)/ (2|(J3;(0))|) 
(blue) and {Jx{t)) / {Jx{0)) (grey). The theoretical fits (full and dashed black line) are based on the parameters d = 55 
(optical density), Fcoi ~ 0.002ms~^ (coUisional rate) and F = 0.193ms^^ (dephasing rate \22\). The rates for driving field 
induced transitions |4,±3) — >■ |4,±4) and |4, ±4) — ^ |4, ±3) are given by fji^T and u^F respectively, where F — 0.002ms~^. 
b) Entanglement ^(t) versus time in ms. Blue data points correspond to the results shown in a). Data points in orange are 
obtained for a lower optical depth (d = 35). The other parameters used in the fits take the same values as in a), ^(t) < 1 
certifies the creation of an inseparable state. The relevant pulse sequence is shown below. The data taken in the absence of 
the driving field (black points) show no entanglement, c) Dissipative entanglement generation in the presence of the pump 
field which incoherently transfers atomic population from undesired levels within F — 4 back to the two level subsystem. The 
pump rate is Fpump = 0.168ms~^. d — 37, the fitting parameters Fcoi and F take the same values as in a) and F — 0.233ms~^. 
The inset shows the evolution of ^{t) after the driving field is switched off after entanglement is generated by dissipation, d) 
Entanglement ^{t) for different initial conditions. The upper curves show a purely dissipative evolution. The lower curves - 
the entanglement generated by dissipation combined with the measurement. Points on the right represent an average over 
measurements of one hour where atoms were kept in a steady state. The used exponential time mode functions are depicted in 
the pulse sequence and are described in the text, e) Schematic illustration of entanglement generation and verification. The 
signal from the detector D for times t > T is used for verification of entanglement in (a-c). In d) the signal taken at f < T is 
given to the verifier as additional information. 



T » For - « 1 and finite 7sT the input- 
output relations reduce to the quantum-nondemolition 
type. The atomic EPR variance at time T can 

be inferred by using the input-output relations ([T]). The 
Stokes operator ^2 (the photon flux difference between 
-1-45° and —45° polarizations with respect to the y-axis 
in Fig. lb) is measured in the time interval [T; T + iprobo] 
(see pulse sequence in Fig. 2b) with the photocurrent 
electronically processed to obtain the relevant light mode: 
yrcadout^ /jT+'p'"'"^ S'°"*(t)cos(m)e-'^=*dt. The param- 
eters K^,7s, (/i — v) of the hnear input-output relations 
are calibrated as described elsewhere [T31[TU]. The atomic 
state reconstruction is calibrated and verified carefully as 
described in detail in the supplemental material SM f3T] , 
where also the modification of the input output equations 
by losses and decoherence is presented. We conclude that 
the measurement of ^ is reliable within the uncertainty 
of ±4% arising from uncertainty in the measurements of 
K^, the detection efficiency 77 and the shot noise of light. 

In the first set of experiments, entanglement is gener- 
ated purely dissipatively. In the first series of this set, 
the pump- and repump fields are turned off at time t = 
(Fig. 2a,b) and the driving (entangling) laser is turned 
on. In the presence of the drive field [P ~ 5.6mW) T2 
is reduced to 6ms and Ti to 34ms. This decoherence 
has been considered the fundamental limitation for the 
entanglement generated by measurements (6j. Here, the 



collective entangling dissipation due to forward scatter- 
ing dominates over the single atom decoherence and leads 
to a rapid reduction of Sj(t) on the time scale of 77^. 
Fig. 2a shows the time evolution of 5]j(t) normalized to 
2|(J^(0))|. For a Coherent Spin State (CSS) ^ = 1, and 
S^css — 2|(Jx)| defines the projection noise (PN) level, 
below which lies the noise level of entangled states. The 
dynamics of 2\{Jx{t))\ due to single atom spontaneous 
emission and collisions on the slow time scale of Ti is 
also shown in Fig. 2a. Fig. 2b presents the time evolu- 
tion of entanglement for two values of the optical depth 
d w 34 (e = 8.5°) and d « 56 (9 = 14.0°). The data is 
well fitted with theory [T^] using the collisional rate Fcoi 
and dephasing rate F j22j compatible with experimental 
values. The details of calculations of the fits are given 
in [21 . The time interval 0.015s over which entangle- 
ment is continuously maintained is several times longer 
than the best previous results obtained for measurement 
induced entanglement O |6] and much longer than T2- 
For comparison, if the driving (entangling) laser is off 
during < i < T and is turned on only at i = T to 
measure the atomic state, ^(T) predictably stays above 
the PN level (black points in Fig. 2b). Also a slight mis- 
match of the Larmor frequencies of the two ensembles 
in the preparation period by ^ 20Hz leads to the disap- 
pearance of the entanglement. This can be viewed as a 
direct consequence of the "which way" information due 
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to the distinguishability of photons emitted by the two 
ensembles. 

In the series presented above, entanglement is created 
in a quasi steady state rather than in a steady state, as 
would be the case for atoms with a true two-level atomic 
ground state, for example in Ytterbium ensembles [53]. 
On the time scale of Ti , atoms are lost to other magnetic 
sublevels of F = 4 and to the level F = 3. This causes 
the eventual extinction of entanglement as described well 
by the theoretical fits shown in Fig. 2a, b with the pump- 
ing rate Fpunip being close to the experimental value. In 
the next series of experiments, the pumping field of an 
optimal strength resonant with the F = A state is kept 
on during the entanglement generation period t > ( 
Fig. 2c). Remarkably, this incoherent process does not 
suppress generation of entanglement, but on the contrary 
brings it further towards a steady state. The entangle- 
ment can now be maintained for 0.04s thanks to pumping 
atoms from sublevels |mi?| < 3 which contribute higher 
noise, back to \mF\ = 4 which is a dark state for the 
pump beam. The eventual loss of entanglement, is in 
part due to atoms which are lost to the F — 3 ground 
state, effectively reducing d. If the entangling mechanism 
is turned off, the entangled state decays in 2ms (inset in 
Fig. 2c), as expected [6j from the decoherence in the dark. 

Finally, we demonstrate generation of steady state 
atomic entanglement. To this end, a repumping field 
F = 3^F = 4is added during the entanglement 
generation, thus closing down the last escape channel 
from the relevant spin system. The atoms reach a steady 
state which is, however, not entangled because the col- 
lective processes are not sufficiently strong to overcome 
the noise added by the incoherent repumping field. The- 
ory predicts [H] that a steady state entanglement can be 
achieved for d = 100, but this is experimentally unfea- 
sible. However, we can use the fact that due to single 
atom decoherence sources, the atomic state is not pure, 
and hence forward scattered light is not completely dis- 
entangled from the atoms. Up to now, measurements on 
light variables have only been used to verify entanglement 
at time T for which only (t > T) have been utilized. 
Using the results of the continuous measurement on the 
open atomic quantum system during the interval t < T, 
we can enhance the entanglement generated by dissipa- 
tion and maintain it in the steady state. In this scenario, 
the verifier (Fig. 2e) receives the classical information 
S'|"*(i < T) which is used to calculate the conditional 
variance var(?/^°"d) = va.r{y'^';f °"'(r) - ayl'ifit < T)). 

Here, y'r'^it < T) = /J' 5^ (0 cos(r!Oe^™*dt 

and the feedback gain a and the time constant 7m are 
optimized to achieve maximal noise reduction. The light 
mode that brings about the best noise reduction is a 
fast growing exponential mode, with 7^ — 0.83ms~^ > 
7 = I/T2 = 0.27ms~^. The conditionally reduced atomic 
variance ^cond = var[i(X/-X//) + i(P/-P//)-^(yf-d + 



2/f°'='^)(t < T))] is then found from var(?/™"d) using the 
same calibrated input-output relations as above. These 
central results are displayed in Fig. 2d, which shows the 
evolution of the variances of the purely dissipatively gen- 
erated atomic state (upper curves) and the entanglement 
produced using the hybrid method including dissipation 
and continuous measurements (lower curves). Each pair 
of curves is taken with two different initial conditions. 
These results demonstrate a very important aspect of 
our work, - they show that the generated steady state is 
independent of the initial state, and that entanglement 
is maintained for up to an hour, if dissipative processes 
are combined with measurements. 

In conclusion, we have observed entanglement of 
macroscopic atomic ensembles generated by dissipation 
and the steady state atomic entanglement. The results 
present a new step in quantum control of entanglement. 
Dissipatively generated entanglement provides not only 
event-ready entangled links for standard protocols but is 
also an elementary resource for future applications in con- 
tinuous quantum information processing schemes, such as 
dissipative distillation and repeater protocols, which al- 
low for the distribution of long-range high-quality steady 
state entanglement [7]. 
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SUPPLEMENTAL MATERIAL 

The first part of the supplemental material contains details 
of the calibration of the atomic spin noise with respect to the 
projection noise and of the orientation of the collective atomic 
spin. 

The second part contains details of the theoretical fits pre- 
sented in Fig.2 of the main text. 



Projection noise calibration 

The variances of the collective atomic operators Pj" -f- P|5 
and X}" — X)° are found from the measured polarization mode 
of the transmitted light yc,3- using the input-output relations 
given in Eq. (1) in the main text, as follows: 

var(P;" -I- Pn)/2 = (var(y°"') - crL(l - k^(^ - i/f))/^'^, 
var(X}" - Xu)/2 = (var(?y°"') - a^„{l - H^i/j. - i/)^))/Kf9.1) 

where (t^„ is the shot noise of light. The normalized EPR 
variance of atomic noise is ^ = var(P|" -I- Pj")/2 + var(X}" — 
XTi)/2. 

Decay of the atomic state with the rate 7oxtra can be in- 
cluded in the input-output equations (Eq. (1)) as follows: 

+eVl - e-^iT ■ Fp+ 



+eVl - e-'^-'T ■ P, 



(S.2) 



with the total decay 7 = l/r2 = 7s +7cxtra, e = 7extr, 
and the coupling constant k = v/(l — e^)(l — e~^T^)/(/i — v). 
The two-cell noise operators (P^ ) = | model a decay towards 
the coherent spin state (CSS), which is a good approxima- 
tion on short interaction timescales. The equations for light 
are adjusted accordingly and used for the reconstruction of 
the atomic noise. The coupling constant is calibrated as 
discussed in [14] . In Fig. |S.l^ , measurements of the coupling 
constant are shown for different atom numbers, monitored by 
the Faraday angle O oc Ja, which is measured by the polar- 
ization rotation of an i ndep endent probe beam propagating 



in x-direction (see Fig. S.2 1, (/i — v) depends on detuning 



and has been found experimentally [19] and derived from the 
theory - both methods give the value 0.16 ± 0.005. 

The imperfect detection efficiency r) = 0.84(3) is included in 
the input-output equations via a simple beam splitter model. 
The light losses in the detection path which are relevant for 
the atomic noise reconstruction have been determined in two 
independent ways. They were measured directly and also in- 
ferred from the atomic noise as a function of the atom number 
presented in Fig. |S.1[ 3. In order to establish the projection 
noise (PN) level we run a series of experiments in which we 
determine ^(t = 0) (more precisely at t = 0.02msec), i.e., 
right after the optical pumping is complete as a function of 
the macroscopic spin Jx- The number of atoms (macroscopic 
spins Jx) is varied by changing the temperature of the cells. 
We use a 1ms drive pulse and perform measurements of 5*2 
as sketched in Fig. |S.2[ We then reconstruct the atomic noise 
with the equations discussed above. 
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Supplemental figure S.l: Calibration measurements: a) 

Normalized initial spin noise variance as a function of O oc 
Jx- The light power was held constant at 5.6mW, the pulse 
duration was 1ms. Red line - linear plus quadratic fit. Black 
line - linear part of the fit presenting the projection noise, b) 
Measurement results for for different Faraday angles 9. 




Supplemental figure S.3; Orientation test for the initial 
atomic state and after 15ms probing with 5.6mW. The time 
evolution of a displacement in the two spin components in the 
rotating frame X oc (Jy) and P oc { Jz) is detected. The black 
dots are the measurement outcomes and the red lines the 
fitted curves, from which the orientation o can be determined. 
The insets show the level structure (with v ~ 20Hz) and the 
pulse sequence. The first probe pulse is optional. 
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Supplemental figure S.2: Measurement setup and geom- 
etry of the experiment: The dc polarization detectors 
measure the Faraday rotation angle B proportional to the 
macroscopic spin Jx- The ac 5*2 detector signal processed by 
the lock-in amplifier (LA) is used to determine the atomic 
quantum spin components Jy,z as described in the main text. 



We also independently measure the degree of spin polar- 
ization, as described below, exceeding 0.992 for all relevant 
atomic num bers. The results of the atomic noise measure- 
StV)) are well fit by ^ oc 9 + 0.004(3)0^ 



ments (Fig 
proves that t 



This 



the quadratic deviation from the PN is small and 
well characterized. For the presented range of number of 
atoms corresponding to = 6 — 11 degrees, ^(0) is a few 
per cent above the PN-level. We also verify that ^ is indepen- 
dent of the duration of the verifying measurement tprobc for 
0.5ms< fprobo <3ms. We conclude that the measurement of 
^ is reliable within the uncertainty of ±4% arising from un- 
certainty in the measurements of k^, the detection efficiency 
?7, and the shot noise of light. 



Measurement of the macroscopic spin orientation 

The optimization of the orientation o = j m-pm, where 
Pm is the population of the magnetic sub-level m, is crucial for 
the success of the entanglement generation. Significant imper- 
fections of the orientation lead to atomic noise above the PN 



level. Also several assumptions in the theoretical model are 
only valid for a highly oriented state. In the experiment the 
orientation was monitored as a function of time during the en- 
tire probe duration using a method closely related to the one 
described in |S2| . After the preparation of the CSS by optical 
pumping, we use a weak RF magnetic pulse at the Larmor 
frequency to create excitations causing a displacement of the 
rotating collective spin. This displacement is read out with a 
weak probe beam. The coherences between the different mag- 
netic sub-levels oscillate at slightly different frequencies due 
to the second order Zeeman shift. In our case this splitting 
is 20Hz. When we look at the evolution of the displacement 
demodulated with an RF frequency, close to the Larmor fre- 
quency SI, a population of sub-levels other than m = 4 will 
manifest itself in a quantum beat signal at multiples of the 
second order Zeeman splitting frequency. In Fig. |S.3[ two 
signals for slightly different orientations are shown and the 
difference is clearly visible. This method is remarkably sensi- 
tive to orientation imperfections. After the pumping, we start 
very close to the CSS with o — 0.998(3). After a 15ms probe 
pulse with 5.6mW the orientation is reduced to o = 0.980(3). 



Comparison of experimental data and theoretical 
predictions 

The generated amount of entanglement ^(t) can be calcu- 
lated as described in detail in [12]. As explained there, the 
produced entanglement is given by 



m = 



S.7(t) + 147V|4.±3)(t) 
N2it){P2{t) + 7) ' 



(S.3) 



where N\ 



|4,±3) 



the number of particles in state 
|4, ±3) (I t)/| ;)) and A^a = iVt + and Pa = | A^t " N^] /N2 
denote the particle number and orientation with respect to 
the two-level subsystem respectively. The time evolution of 
the EPR variance Sj = var {Jyj — Jy.ii) + var {Jz.i ~ Jz,ii) 
(see main text) can be calculated using Eq. (21) in [12] if the 
collective decay rate (dP) and the effective dephasing rate 
(P) as well as N2 (t) and P2 (t) are known. In the following we 
describe how these quantities can be inferred and explain the 
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theoretical fits to the measured data presented in the main 
text. 

The essential features of the experiment can be described by 
means of a simplified model, which allows one to take addi- 
tional dynamics due to the multilevel character of Cesium 
into account and involves only the three atomic states |4, ±4), 
|4, ±3) and = |3, ±3). For the timescales considered 

here, atomic population in other sublevels can be neglected. 
Since the model is primarily intended to describe qualita- 
tively the physical effects observed in the experiment with 



very few parameters, we use r|4_-|-4)_ 



and r 



■\h) 



|4,±3)-!.|h> 



|h>^|4,±4) - 

d ~ , . 



"|h>^|4,±3) = Tin such that 

-(rout + 2rin)7V2W + 27vri, 
-(r3,4 + r4,3 + rout)P2W 

+ (r3,4-r4,3)iV2(t)/7V, 



where P2 (t) = P2 {t)N2 {t) /N . Here and in the following we use 



the abbreviations F 



4,±4>->|4,±3) 



r4,3 and Fm^ 



±3)^|3,±4> 



Fa, 4. Atomic transitions can be either induced by the driving 
field or due to collisions. Since the thermal energy of atoms 
is much larger than the atomic level splittings, we assume the 
same coUisional rate Fcoi for all atomic transitions. Accord- 
ingly, 



L3,4 
^4,3 



= F + Fcoi, 



out — i L + J- col, 
Fin ~ Fcol, 



where the abbreviations F 



|4,±4)^.|4,±3) 



|4,±3>->|3,±4> 



3,4 



have been used. 



— F4,3 and 
and v^V are the 



driving field induced cooling and heating rate respectively. 
Fl"* is the rate at which atoms leave the two level subsystem 
due to radiative transitions caused by the driving field. The 
number of free parameters in these equations can be reduced 
to two, using the experimentally determined time derivative 
of the atomic polarization P — {Jx{t)} / {Jx{0)) at time t — 



Pit) 



-(F4.3+4Fout)A^|4,±4)(0) + (F3.4-3Fout)iV|4,±3)(0) 



where it is taken into account that the initial spin state is not 
perfectly polarized, but contains a small fraction of atoms 



in state |4, ±3). The initial populations A'"|4 ±4) (0) = 0.99, 
-^|4,±3>(0) = 0.01 and Ai'h(O) = are estimated based on 
measurements of the orientation of the initial spin state af- 
ter optical pumping which are described in the first part of 
SM. Using this constraint, P(t) can be fitted with two free 
parameters. This way, fixed expressions for P2{t) and N2{t) 
are obtained. As mentioned above, the values of the collec- 
tive decay rate dF and the dephasing rate due to noise effects 
F jS3j have to be known in order to calculate the generated 
amount of entanglement as described in [12]. These param- 
eters are determined from the experimentally obtained slope 
of the variance Ej(t)/ (2|(Ja;)|) at time t = 



d Sj(t) 
dt2\{J40))\ 



= -4iVdFP2(0)(l-P2(0)/(M-/^)') 



+ 7FinP2 (0) - 7 (Fout +F3,4 -F4,3)iV|4,±3> (0), 



and the decay of the transverse spin {Jy{t)} = 
g-i(f+drP2(t))t^j^(0)), where Eq. (21) in [12] and the 
identities relating quantities defined with respect to a 
two-level system to quantities defined with respect to a 
multi-level structure presented in Sec. IV. A in [12] have been 
used. 



[SI] In the actual experiment, the 7cxtra will collect all decay 
mechanisms coming from spontaneous emission, collisions 
and magnetic field inhomogeneity 

[S2] B. Julsgaard, J. Sherson, J. S0rensen, and E. S. Polzik, 
J. Opt. B, 6, 5 (2004). 

[S3] In the absence of pump or repump fields, the dephasing 
rate associated with noise effects F = F3,4 + F4,3 -I- Fdeph 
includes the effective single particle cooling rate F3.4 = 
/i^F -I- Fcoi, which is a sum of the rates corresponding to 
cooling processes due to the driving field and collisions, 
the effective single particle heating rate F4,3 = z/^F -I- Fcoi 
as well as pure dephasing Fdeph- 



